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ABSTRACT 

Recently gathered observational data on a sample of Type la Su- 
pernovae (SNe la) reveal a wide distribution of expansion velocities of 
the Fe cores, measured from the width of the nebular lines. Moreover, 
the velocity appears to correlate with the luminosity decline rate af- 
ter maximum light, Ami5(5). Since it has been shown that for SNe la 
Ami5(5) correlates with the absolute magnitude at maximum, this then 
implies a relation between the expansion velocity of the Fe nebula and 
the luminosity at maximum. 

Physically, the maximum luminosity is related to the mass of syn- 
thesized ^^Ni, whereas the FWHM of the lines is related to the kinetic 
energy of the ejecta. Our finding constitutes observational proof of the 
theoretical prediction that the two quantities have to be related. 

Subject headings: supernovae: general 
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1. Introduction 

Type la Supernovae (SNe la) are thought 
to result from the explosion of mass-accreting 
white dwarfs (WD). Since SNe la are bright, 
they can be used to determine distances if 
their absolute brightness can be calibrated. 
The fact that all SNe la seemed to show sim- 
ilar light curves and spectra encouraged their 
use as standard candles. In addition, SNe la 
are important for galaxy chemical evolution 
since they provide most of the Fe in the uni- 
verse. 

With the accumulation of new and more 
accurate data, the paradigm that SNe la are 
all similar has been questioned. In particu- 
lar, the discovery of objects such as the very 
bright SN 1991T (Filippcnko ct al. 1992a; 
Phillips et al. 1992; Ruiz-Lapuente et al. 
1992) and the very faint SN 1991bg (Filip- 
penko et al. 1992b; Leibundgut et al. 1993; 
Turatto et al. 1997) showed that SNe la ex- 
ist with very different properties from those 
of the classical objects defining this class, e.g. 
SNe 1981B (Branch et al. 1983) or 1992A 
(Kirshner et al. 1993). Other peculiar SNc la, 
both over- and under-luminous, have subse- 
quently been discovered, showing that such 
events, although less frequent than 'normal' 
SNe la, are certainly not just rare excep- 
tions. The progenitor masses and the ex- 
plosion mechanisms leading to these peculiar 
SNe la are currently topics of debate, but it 
is clear that the over-luminous objects pro- 
duce more ^^Ni than 'normal' SNe la, while 
the under-luminous ones produce less (Maz- 
zah et al. 1993, 1995, 1997). 

Objects such as SNe 1991bg and 1991T are 
not only different in luminosity from the aver- 
age SN la, but their spectra are also different, 
at least at phases before and around maxi- 



mum (however, Phillips et al. (1992) have 

shown that SN 1991T 14 days after maximum 
was indistinguishable from other SNe la), re- 
fiecting differences in temperature and also in 
the composition of the outer layers. It has 
been argued that by selecting SNe la on the 
basis of a 'normal' spectrum at about max- 
imum or earlier, one can define a subgroup 
of 'normal' SNe la which can be used as stan- 
dard candles (Branch & Miller 1993; Vaughan 
et al. 1995; Saha et al. 1997). The vast 
majority of well-observed SNe la belong to 
this group, and so peculiar objects such as 
SNe 1991bg and 1991T can be safely removed 
from the sample. 

Therefore, if the 'spectroscopically normal' 
SNe la are to be used as standard candles, all 
of their properties should be similar. How- 
ever, observational evidence has been accu- 
mulating for the existence of a relation be- 
tween the absolute magnitude at maximum 
of SNe la (normal and peculiar) and the 
rate at which their light curves decline af- 
ter maximum. This is known as the Phillips 
or M — Amis relation, which indicates that 
the SNe whose rate of decline is slower are 
those that reach a brighter absolute magni- 
tude (Phillips 1993). Indeed, the Hubble di- 
agram for a sample of distant SNc la has a 
much smaller dispersion when the M — Amis 
effect is corrected for (Hamuy et al. 1996b). 

Model hght curves have shown that the ori- 
gin of the Mb — Amis relation could lie in the 
existence of a range in the Ni mass and dis- 
tribution produced in the explosion (Hofiich 
et al. 1996; Cappellaro et al. 1997). The 
amount of Ni produced is directly related to 
the energy of the explosion and is reflected 
in the kinetic energy of the ejecta. Thus we 
address the question: can any direct observa- 
tional evidence for this scenario be found? 
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2. The late-time spectra 

We have analyzed the photometric and 
spectroscopic data on SNe la collected at ESO 
over the last several years. We concentrated 
on spectra obtained at epochs of about one 
year after maximum. At epochs greater then 
100-150 d, SNe la display a nebular spectrum: 
the ejecta have become so thin that they are 
transparent to radiation. The gas is heated 
by the deposition of some of the energy car- 
ried by the 7-rays and the positrons produced 
by the decay of ^^Co (the daughter nucleus 
of ^^Ni) into stable ^^Fe, while cooling takes 
place via forbidden line emission. Forbidden 
lines of Fe II and Fe III dominate the spec- 
trum, confirming that iron-group elements are 
prevalent in the inner part of the ejecta. The 
outer layers, where lighter elements are dis- 
tributed, have too small a column density to 
absorb the surviving 7-rays (the positrons, 
which actually dominate the energy input af- 
ter about 200 days, deposit essentially in situ) 
and thus do not contribute significantly to the 
spectrum. 

A representative sample of SNe la spectra 
at an epoch of about 1 year taken from the 
ESO+Asiago SN archive is shown in Fig.l. 
Although the spectra shown in Fig.l belong 
to SNe that differ by up to 3 magnitudes 
in brightness at maximum (from SN 1991bg 
to SN 1991T), it is interesting to note that 
their nebular spectra are quite similar, ex- 
cept for the line widths. In particular, the 
two strongest lines in the blue have similar ra- 
tios, when reddening has been corrected for. 
Since the bluer feature (4700A) is almost ex- 
clusively due to Fe III forbidden lines, while 
the red feature (5300A) is formed by a blend 
of Fe II and Fe III forbidden fines, their ratio 
is sensitive to the ionization/excitation condi- 
tions (and hence temperature and density in 



the nebula). Preliminary modelling work on 

the nebular spectra (Mazzali et al. 1998, in 
prep.) indicates that the observed range of a 
factor 1.5 in the line ratio can be obtained 
for a difference in temperature and ioniza- 
tion of not more than 10%. This must im- 
ply that the energy deposition rate per unit 
volume in the SN nebula is quite similar in 
all SNe. Indeed Cappellaro et al. (1997) 
showed that the difference in brightness at 
maximum is preserved, or even increased at 
the later epochs. Since our sample includes 
objects from the least luminous SN la known 
to date, SN 199 Ibg, to the most luminous one, 
SN 1991T, the entire range of SNe la appears 
to share these same properties. 

The composition of the emitting nebula ap- 
pears to be similar in all SNe: all spectra are 
dominated by lines of Fe and Co, although the 
mass of these elements produced is probably 
very different in the various objects, ranging 
from about 0.1 in SN 1991bg (Mazzah et 
al. 1997) to approximately 0.6 M0 in 'normal' 
SNe la and about 1 Mq in SN 199 IT (Spy- 
romilio et al. 1993). This confirms that only 
the inner Fe sphere is significantly affected by 
the heat input following the absorption of the 
7-rays and the positrons, essentially a local 
process, and that mixing is not significant, at 
least at the lowest velocities. 

The second distinguishing property of the 
nebular spectra of Fig.l is that the fines have 
different widths in different SNe. The line 
widths span about a factor of 5: the faint 
SN 199 Ibg has the narrowest lines, and the 
bright SN 1991T has the broadest lines. 

If the physical conditions (temperature, 
density) in the nebulae are similar, as indi- 
cated by the line ratios, and if the brighter 
SNe produce more ^^Ni, it should be ex- 
pected that the more massive Ni sphere in 
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the brighter and more energetic SNe is larger 
than in the fainter SNe. Since the expansion 
of the ejecta is homologous (i.e. f oc r), a 
correlation between the width of the nebular 
lines and the brightness of the SN should be 
expected. 

Excluding the case of SN 1991bg, where 
several individual lines are resolved, in all 
other cases the observed lines are the result of 
very strong blending. In particular, the fea- 
tures in the range 4500-5500A can be fitted 
with three gaussians to measure position, in- 
tensity and FWHM (the continuum being set 
to zero). The FWHM of the strongest feature, 
measured at ~ 4700A, was then selected as 
the reference value, because it is always nar- 
rower than that at ~ 5300A and has a much 
better S/N than the much weaker feature at 
~ 5000A. 

The value of FWHM is representative of 
the ejecta velocity, but a knowledge of the 
density /velocity distribution would be needed 
to translate it into a kinetic energy. In ad- 
dition, when blending is severe, the value of 
FWHM is an overestimate of the velocity of 
the Fe nebula. 

The effect of line blending can be esti- 
mated using synthetic spectra (Mazzali et al. 
1998, in prep.). By comparing the velocity 
used in the model for the Fe-ncbula with the 
FWHM of the ~ 4700A feature measured on 
the model synthetic spectra we found that 
blending can be neglected only for input ve- 
locities of up to 3000 km/s. The relation be- 
tween input velocity and measured FWHM 
has a discontinuity when blending sets in, be- 
tween about 3000 and 6000 km/s, but for 
higher velocities a monotonic relation is re- 
covered. The velocity of the model Fe-nebula 
is marked in the upper x-axis of Fig. 2. 

The measured FWHM should be plotted 



against the SN absolute magnitude. Unfor- 
tunately, the absolute magnitude is a very 
uncertain quantity, since it depends on the 
adopted distance scale and extinction correc- 
tion. To avoid these problems, we make use of 
the relation between the absolute magnitude 
at maximum and Amis (5), the difference in 
B magnitude between the time of maximum 
and 15 days later (Phillips 1993, Hamuy et 
al. 1996b). This latter quantity can be mea- 
sured directly, so that two observables can 
be compared. Measurements of Amir,{B) are 
readily available in the literature for most of 
the objects in our sample. For SNe 1994ae 
and 1995D, we measured Ami5{B) on the 
hght curves pubhshed by Riess et al. (1996), 
whereas for SNe 1991M and 1993L we used 
unpublished ESO photometry, which unfor- 
tunately does not cover the time of maximum 
very well. We also note that the hght curve 
of SN 1989M is very poorly sampled. 

A plot of FWHM(4700A) v. Arm^iB) is 
shown in Figure 2. It is apparent that a cor- 
relation exists between the two quantities (the 
Spearman rank-order correlation coefficent is 
—0.81). Note also that a significant spread 
in both line width and decline time is present 
even among 'normal' SNe la. Even when the 
extreme case of SNe 1991bg is removed from 
the sample, the null hypothesis of no correla- 
tion between Ami^iB) and FWHM(4700A) 
has a probabihty of only 0.003. Furthermore, 
the FWHM measured for SN 1991bg is lower 
than for other SNe both because of the intrin- 
sically low velocity and the absence of blend- 
ing. 

3. Discussion 

Our results establish a connection between 
the luminosity of the SNe and their kinetic 
energy. Since the maximum luminosity of 
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a SN la depends essentially on the mass of 
^^Ni produced, the observed difference in line 
width and strength in the late-time spectra 
is a strong indication that SNe la produce a 
range of Ni masses. In the case of the peculiar 
SNe la, the properties at early times are suffi- 
ciently different that the spectra are affected, 
while for 'spectroscopically normal' SNe this 
is not the case, and the differences emerge 
only in the rate of decline and in the nebular 
line width. However, Branch &: van den Bergh 
(1993) noted a sequence in the velocity of the 
Si II lines in SNe la about 10 days after max- 
imum. This relation to galaxy type may be 
a less quantitative manifestation of the same 
phenomenon we are witnessing at late time. 
Also, Nugent et al. (1995) suggested a 'spec- 
troscopic sequence' among SNe la near max- 
imum. All the SNe in their sample are also 
included in ours, and the ordering of the ob- 
jects by nebular line width is similar to what 
they suggest based on the value of V{Max). 
They also hypothesized that the sequence is 
the result of the different amounts of ^^Ni pro- 
duced by the various SNe. 

We can assume that the line width de- 
fines the volume of the Ni sphere, and that 
this translates directly into a Ni mass. If we 
accept the estimate of the Ni mass for SN 
1991T (- 1.0 Mq, Spyromilio et al. 1993), 
we can scale the other SNe accordingly. We 
find then that SN 1986G should have a Ni 
mass of ~ 0.3 M©; objects such as SNe 1981B 
and 1992A, which we could call 'faint nor- 
mal', have ~ 0.4 — 0.5 Mq, while 'bright 
normal' objects like SN 1990N and 1994D 
have ~ 0.5 — 0.6 M©, and appear to be the 
only SNe producing the 'canonical' amount 
of ^^Ni. These values are only rough esti- 
mates, and further modelling work is required 
to obtain more accurate results. Note that 



Ruiz-Lapuente & Lucy (1992) determined for 
SN 1986G Mi^^m) = 0.38 ± 0.03 Mq, but 
they used a larger reddening. Finally, the 
very subluminous SN 199 Ibg appears to have 
a very small Ni mass, of about 0.01 Mq. This 
is smaller than the pubhshed values (Mazzali 
et al. 1997), but note that in the case of 
SN 1991bg blending does not affect the line 
width, and that the ratio of the fluxes in the 
blue and the red feature is the smallest among 
all the SNe in our sample, indicating a lower 
degree of ionization and therefore a smaller 
fractional abundance of the Fe III ion, which 
is responsible for the blue feature. The mod- 
els for SN 1991bg (Mazzah et al. 1997) indeed 
require a higher Ni density than in other SNe. 
This leads to a lower ionization degree, and 
the flux is emitted more in the red line near 
5300A. 

Models for the explosion of WD's lead- 
ing to SNe la which produce less than the 
canonical 0.6 Mq of ^^Ni typically require pro- 
genitors less massive than the Chandrasekhar 
limit of 1.4 Mq (see e.g. Livne & Arnett 1995 
and references therein). Hoflich et al. (1995) 
presented Chandrasekhar mass, pulsating de- 
layed detonation models for SNe la. They can 
vary the mass of ^^Ni produced by varying 
the density at which the transition from de- 
flagration to detonation occurs. Models that 
synthesize less ^^Ni produce more Si-group el- 
ements, so that all models produce about the 
same kinetic energy. They claim this family 
of models can explain the observed differences 
among SNe la light curves. Our results prove 
that the kinetic energy is actually related to 
the ^^Ni mass, so that burning to NSE is nec- 
essary in all realistic explosion models. Since 
the rapid decline of the faint SNe la does not 
appear to be the result of a large kinetic en- 
ergy in a Chandrasekhar-mass explosion, it 
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is then likely to be due to a small exploding 
mass. The fact that SN 1991bg was the re- 
sult of the explosion of a sub-Chandrasekhar 
WD has received support from several inde- 
pendent studies. If the range of ^^Ni masses 
which we observe in the nebular spectra is 
indeed the result of a range of progenitor 
masses, our understanding of the physics and 
the nature of SNe la would be revolutionized. 

Formally, we can fit the data of Fig. 2 with 
a straight line. Since errors affect both axes, 
we performed a weighted, linear least square 
fit (Press et al. 1992), which gives the line 
shown in Fig. 2. However, the variance of the 
points around the line seems too large to be 
due only to the estimated errors (the reduced 

= 4.1, and even excluding SN 1991bg we 
get = 2.54) 

Whereas a possible interpretation of the 
large variance of the points in Fig. 2 is that 
the relation Ami5(5) vs FWHM^oQd is not 
linear, it is possible that we are seeing an in- 
dication of the existence of a 'second parame- 
ter' determining the shape of the light curves 
(or possibly of the velocities). 

A similar suggestion was recently made by 
Hoflich et al. (1996) in the context of a discus- 
sion of the My — Am2o{V) relation. As the 
same authors commented, however, the un- 
certainties in the distance moduli and in the 
extinction corrections, and hence in the cal- 
ibration of the SN absolute magnitude, cast 
some doubts on the reality of their finding. 
The variables plotted in Fig. 2 are free from 
such uncertainties, and therefore the result is 
particularly intriguing. This 'second param- 
eter' could be related to the mode of explo- 
sion, and possibly also to the difference in the 
masses of the progenitors through the ratio of 
the mass of synthesized ^^Ni and the mass of 
the ejecta (Cappellaro et al. 1997). The ex- 



istence of a 'second parameter' might weaken 
the potential of using the decline rate as a 
means to determine SN luminosities. 

However, the nebular spectra may offer an- 
other possibility. From a theoretical point of 
view, the absolute magnitude at maximum 
should be directly correlated to the mass of 
synthesized ^^Ni, with an associated uncer- 
tainty depending on the explosion model. If 
the ®^Ni mass can be obtained from the neb- 
ular spectra, we may be able to use the late- 
time spectra of SNe la as distance indicators, 
thus avoiding the uncertainties introduced by 
the 'second parameter'. If a relation can be 
established between the width of the nebu- 
lar lines and the absolute magnitude at max- 
imum, and once it is well calibrated with the 
acquisition of more data, the absolute mag- 
nitude of a SN, and hence its distance, could 
be determined straightforwardly from a line 
width measurement on a nebular spectrum at 
an epoch of about one year. Precise timing of 
the spectra is not so important, since we do 
not expect and have not seen any significant 
systematic evolution of the line width while 
the SNe are in their nebular phase. Also, the- 
oretical models can be used to derive the ^^Ni 
mass (Ruiz-Lapuente & Lucy 1992) and hence 
the luminosity at maximum (Cappellaro et al. 
1997). Against this however is the fact that 
the SNe at one year are at least 7 magnitudes 
fainter than at maximum light. 

If SNe la produce a continuous range of 
^^Ni masses and kinetic energies, it is also 
quite likely that the masses of the progenitors 
vary, thus favouring the sub-Chandrasekhar 
scenario for a large number of SNe la. Even 
objects that have been regarded as normal be- 
cause of their spectroscopic appearance may 
indeed be subluminous and submassive, al- 
though less so than the well known cases. 
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which are only the extremes of a continuous 
distribution. 

This work was conducted in the framework 
of the ESO Key Programme: Supernovae. It 
is a pleasure to thank L.B.Lucy for several 
useful duscussions and comments on an ear- 
lier version of this manuscript. 
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Fig. 1. — Optical spectra in the region 4400- 
5700A for four representative SNe la. The 
spectrum of SN 1986G has been corrected for 
extinction assuming E(B-V)=0.6. The dotted 
lines show the multiple gaussian fit performed 
to measure the FWHM of the emission. For 
SN 1991bg, the dashed line is the gaussian fit 
of the 4660A hne. 

Fig. 2. — The value of Ami5(S) is plotted 
against the FWHM of the ~ 4700A emission 
feature. The best fit line is a weighted, least 
squares fit of the data points, and has equa- 
tion Ami5(S) = (1.24±0.03)-(0.10±0.01)x 
{FWHM -15000). 
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Table 10 
The sample of SNe Ia 



SN 


epoch 
days 


^max 


AmB(15) 


Av 


ref 




ref (meth) 


Mv 


FWHM 
(10^ km s-i) 


ref 


1981B 


300 


11.93 


1.10 ±0.05 


0.00 


a 


31.10 


1 (Ceph) 


-19.17 ±0.15 


14.1 ±0.5 


n 


1986G 


300 


11.38 


1.73 ±0.07 


1.95 


a 


27.86 


a (SBF/PNLF) 


-18.43 ±0.32 


12.2 ±0.4 


o 


1989B 


330 


11.95 


1.31 ±0.07 


1.05 


b 


30.28 


1 (Ceph) 


-19.38 ±0.18 


12.6 ± 1.9 


P 


1989M 


300 


12.30 


1.10 ±0.20 


0.11* 


c 


31.10 


c (T-F) 


-18.91 ±0.36 


17.7 ±0.4 


E 


1990N 


245 


12.70 


1.03 ± 0.06 


0.00 


d,a 


32.03 


I (Ceph) 


-19.33 ± 0.23 


15.4 ±0.8 


q 


1991M 


150 


14.10 


1.50 ±0.10 


0.05* 


E 


32.72 


# (v3K) 


-18.67 ±0.60 


14.0 ±0.7 


E 


1991T 


282 


11.51 


0.95 ±0.05 


0.39 


d,m 


30.60 


b (T-F) 


-19.48 ±0.30 


18.0 ±0.5 


m 


1991bg 


203 


13.96 


1.95 ±0.10 


0.15 


f 


31.65 


i (GCLF) 


-17.84 ±0.36 


3.5 ±1.0 


f 


1992A 


300 


12.55 


1.47 ±0.05 


0.00 


a 


31.35 


i (GCLF) 


-18.79 ±0.16 


14.8 ±0.4 


E 


1993L 


269 


13.20 


1.50 ±0.10 


0.03* 


E 


31.64 


# (v3K) 


-18.47 ±0.60 


13.3 ±0.6 


E 


1994D 


300 


11.90 


1.26 ±0.05 


0.09 


g,h 


30.86 


a (SBF) 


-19.05 ±0.12 


15.3 ±0.4 


E 


1994ae 


182 


13.05 


1.00 ±0.10 


0.42 


e 


31.28 


c (T-F) 


-18.65 ± 0.41 


14.6 ±0.6 


E 


1995D 


363 


13.40 


1.00 ±0.10 


0.21 


e 


32.73 


# (v3K) 


-19.54 ±0.40 


15.2 ± 1.0 


E 


1996X 


299 


13.25 


1.30 ±0.10 


0.20* 


E 


32.77 


# (v3K) 


-19.73 ±0.40 


15.0 ±0.7 


E 



References. — E = ESO+Asiago SN arcliivc; a = Hamuy ct al. 1996a; b = Phillips 1993; c = Vaughan ct al. 1995; d = Lira et al. 
1998; e = Riess et al. 1996; f = Turatto ct al. 1996; g = Patat ct al. 1996; h = Ho & Filippenko 1995; i = Delia Valle et al. 1998; 1 
= Saha ct al. 1997; m = Mazzali ct al. 1995; n = Branch 1984; o = Cristiani ct al. 1992; p = Wells et al. 1994; q = B. Leibundgut 
1996, priv. comm.; # calculated from v3k from RC3 adopting H=65 km Mpc~^ 

*only galactic absorption from RC3 
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